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ABSTRACT 

An experimental wind-tunnel study to 
determine the divergence characteristics of a 
high-aspect ratio, forward-swept wing has been 
conducted in the NASA Langley Research Center 
(LaRC) Transonic Dynamics Tunnel (TOT). The 

rectangular wing used for this study had a panel 
aspect ratio of 9.16 (A = 0°) and the sweep 
angle could be set at A = 0°, -15°, -30°, -45°, 
or -60°. A rectangular winq tip shape was 
tested at each of these sweep angles. In 

addition, a tip shape parallel to the freestream 
flow was tested for a wing sweep angle of 
A = -45°. The root of the wing was cantilever 
mounted to the wall of the wind tunnel. 
Divergence conditions were measured at M = 0.4 
for each sweep angle and tip configuration 
tested Subcritical response techniques were 
used to extrapolate to the divergence conditions 
during the wind-tunnel test. The primary 
objective of this test was to obtain data which 
could be used to verify for this configuration 
the divergence prediction capability of an 
aeroelastic analysis code. Subsonic lifting 
surface theory (kernel function) aerodynamics 
are utilized by this particular code. The 
analytical predictions of divergence were found 
to be significantly conservative at all forward 
sweep angles. At A = -45°, the analysis was 14 
percent conservative. The effect of the two tip 
shapes on the divergence dynamic pressure was 
predicted accurately by the analysis. The 
divergence condition for the tip shape parallel 
to the flow occurred at a dynamic pressure 14 
percent higher than the divergence condition 
with a rectangular tip shape. 


NOMENCLATURE 

AR aspect ratio 

El bendinq stiffness, lb- ft 2 

f frequency, Hz 

g incremental damping 

GJ torsional stiffness, lb* ft 2 

M Mach number 

M e moment about elastic axis, ft* lb 

q dynamic pressure, lb/ft 2 

% diveraence dynamic oressure predicted 

by a subcritical response technique, 

ib/n 2 

^ref reference dynamic pressure, Ib/ft 

V velocity, ft/sec 

y spanwise distance along elastic axis 

from wing root, ft 
a root angle of attack, deg 

a divergence index parameter = l-fa/q^) 

X slope of M e -versus-a curve, ft*lb/dea 

A md-chord sweep angle, deq 

p density, lb*sec 2 /ft 4 


Subscripts 

a analysis result 

n nth value 

r rth value 

45 forward sweep angle, deg 


INTRODUCTION 

The combination of vertical take-off and 
high-speed flight is a desirable feature for 
many aircraft applications. The X-Wing concept 1 
vehicle, currently under examination in a joint 
NASA/DARPA flight test demonstration program, is 
an excellent example of a vehicle which has both 
of these capabilities. The X-Wing achieves 
these features by flying in either a helicopter 
or a fixed-wing mode. The "conversion" from 
one configuration to the other takes place 
during flight by stopping/starting the four- 
bladed rotor. In the fixed wing configuration, 
the blades (or wings in this instance) are 
oriented at A = +45° and A = -45°. Therefore, 
two of the rotor blades are actually high-aspect 
ratio, forward- swept wings. The scope of the 
current flight test program is to demonstrate 
the conversion phase of flight. The Rotor 
System Research Aircraft (RSRA) 2 will be 
utilized to achieve this first step toward a 
successful X-Wing vehicle. 

The present study was conducted to obtain 
data which could be used to verify the 
prediction capabilities of a currently available 
aeroelastic analysis code for this unusual 
configuration (high-aspect ratio, forward-swept 
wing). An experiment was conducted in the NASA 
Langley Transonic Dynamics Tunnel (TDT) on a 
previously existing model with an aspect ratio 
similar to that of the RSRA/X-Wing rotor blade. 
This model had a panel aspect ratio of 9.16 
(unswept) compared to the RSRA/X-Wing design 
value of 9.61 (see Fig. 1). Aeroelastic 
analyses were conducted for each condition 
tested in the TDT for this comparison The 
wind-tunnel model was tested at various forward 
sweep angles to determine the effect of other 
azimuth angles which rotor blades on a X-Wing 
vehicle would encounter during the conversion 
phase of flight. A rectangular tip shape (as 
currently planned for the RSRA/X-Wing vehicle) 
was used during most of the experiment. A tip 
parallel to the flow in the A = -45° position 
was also tested for further correlation with 
analysis. General aeroelastic characteristics 
of the high-aspect ratio, highly swept wing and 
the prediction capabilities of the analysis code 
are discussed in this report. 


TEST APPARATUS AND PROCEDURES 


Wind-Tunnel Model 

The model used in this study was untapered, 
had a 4.51 ft semi-span, and had a semi-span 
aspect ratio of 9.16 in the unswept position as 
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Fig, 1.- Comparison of planform Geometry of the 
X-Wing flight configuration and the 
wind-tunnel model. 


shown in Fig. 1. The airfoil section was a NACA 
0014. The model wing was constructed of a 
layered fiberglass shell, which provided both 
structural stiffness and the airfoil shape, with 
a rectangular aluminum spar located at the 30 
percent chord position to increase the bending 
stiffness. This construction, and the airfoil 
shape, are illustrated in Fig. 2. Semicircular 
wing tips made of balsa wood were used to 
improve the flow over what otherwise would have 
been a blunt wing tip in the forward-sweep 
posi tions. 

A splitter plate was mounted to the support 
system at the wing root to provide a reflection 
plane for the model. Fig. 3 is a photograph of 
the model mounted in the tunnel in the A = -45° 
position. The model was clamped in a canti- 
levered manner to the support system which was 
clamped to the wind-tunnel sidewall turntable. 
This allowed the angle of attack to be changed 
during testing. 

The wing was positioned manually to sweep 
angles of 0°, -15°, -30°, -45°, or -60°. A 
composite photograph showing the model in the 
various sweep positions is shown in Fig. 4. The 
model was tested with a rectangular wing tip at 
each azimuth angle as shown in the figure. In 
addition, the wing tip was modified such that 
the tip was parallel to the freestream flow when 


Two layers of 0.004" thick 
fiberglass skin 

NACA 0014 airfoil 
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0.375" steel bail spaced 
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for mass distribution 



Fig. 4.- Composite photoaraoh showing each sweep 
angle tested. 


tested in the A = -45° position. Sketches of 
the two tip shapes are shown in Fig. 5. 

Vibration and stiffness characteristics of 
the wind-tunnel model were measured. Natural 
frequencies were measured from the output of 
electrical resistance strain-gage bridges which 
were mounted near the wing root and oriented 
such that bending and torsional strains were 
measured. Measured and calculated natural 
vibration frequencies are given in Table 1. 


Wind Tunnel 

The model tests were conducted in the 
Langley TDT 3 . The TDT is a continuous-flow wind 
tunnel capable of subsonic and transonic 
speeds. The test section of the tunnel is 16.0 
ft square with cropped corners. The TDT can use 
either air or Freon* 12 as the test medium. 


Fig. 2.- Cross-section view of model 
construction. 


*Freon: Registered trademark of E. I. duPont de 

Nemours & Co., Inc. 
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Table 1* Natural vibration frequencies of the 
wind-tunnel model. 


Mode 

Calculated 
f, Hz 

Measured 
f, Hz 

1-First bending 

5.4 

5.3 

2 -First chordwise hendina 

22.0 

23.8 

3-Serond hendinq 

33.6 

33. n 

4-First torsion 

62.5 

71.3 

5-Third hendina 

92.3 

89.5 



Fiq. 5 - Tio shapes tested in the A = -45° 
configuration. 


Freon 12 was used for the present study. A Mach 
number range from near zero to 1.2 and a 
pressure range from near zero to atmospheric 
pressure are available. 


Test Procedures 

Experimental predictions of static 
aeroelastic divergence were made using 
subcritical response techniques 4 . For each 
sweep anale tested, subcritical data were taken 
at gradually increasing values of dynamic 
pressure At each dynamic pressure the model 
angle of attack was first adjusted to a 1-g lift 
condition so that the weight of the model was 
aerodynamical ly supported. The angle of attack 
was then incrementally increased and the root 
bending moment was measured at each angle of 
attack. These data were used to predict the 
dynamic pressure at which divergence would 
occur. 

Two subcritical response techniques, an 
improved static Southwell method and the 
divergence index method, were used during this 
test. Reference 4 discusses these prediction 
methods in greater detail. The improved static 
Southwell method uses all of the angle of attack 
data obtained at each dynamic pressure tested 
(Fig. 6). The slope, X, of the M e -versus-a 
data is plotted against X/g for each dynamic 


pressure (Fiq. 7). The slope of a linear fit of 
the resulting plot is the predicted divergence 
dynamic pressure. 

The divergence index method also uses all 
of the angle of attack data obtained at each 
dynamic pressure to calculate X. The divergence 
index parameter, 


A = 1 - (q /q ) 
n 'W, 

1 ~<X n /xJ 
n r 

is then calculated. The subscript r indicates a 
reference condition, which for this test was the 
lowest dynamic pressure data set obtained at 
each sweep angle. Using the divergence index 
method, A is plotted versus q (Fig 8) and a 
linear fit is imposed on the experimental 
results such that A is related to q as follows 

A = 1 - q 
9D 

Therefore, at a dynamic pressure of zero, the 



(a) A = -45°, M = 0.4 



(b) A = -60°, M = 0.4 

Fig 6.- Typical data obtained for subcritical 
response divergence predictions. 
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(h) A = -60°, M = 0.4 

Fio. 7.- Tymcal results of the Southwell method 
divergence prediction technique. 


diveroence index parameter eouals one Further, 
at A = 0 the straight line extrapolation 
intercepts the predicted divergence dynamic 
pressure. 


ANALYTICAL TOOLS 

Aeroelastic analyses were conducted for the 
wind-tunnel model to determine the validity of 
the analysis code for a hiqh-aspect ratio, 
highly-swept wing. A vibration analysis was 
performed with the finite element method 
program 5 Engineering Analysis Language (EAL). 
General beam elements were utilized to assemble 
the finite element model. The beam element 
simulation for the EAL analysis is shown in 
Fig. 9. Elements were arranged to simulate the 
stiffness properties along the measured elastic 
axis of the model. Constant stiffness 
properties derived from laboratory measurements 
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q, lb/ft 2 


(b) A = -60°, M = 0.4 

Fiq. P.- Typical results of the diveroence index 
method prediction technique. 


were used for the beam elements. The measured 
stiffness properties and the constant property 
values used in the finite element model are 
shown in Fig. 10. Lumped masses were located to 
simulate the chordwise center of gravity 
position of the model. These lumped masses were 
offset from the structural elastic axis by very 
stiff bean elements to virtually eliminate any 
chordwise warpage of the analytical model. 
Additional elements (with no mass) were located 
at the leading and trailing edges to obtain 
displacements around the perimeter of the model 
for use in the aeroelastic analysis. Natural 
frequencies, natural mode shapes, and 
generalized masses were calculated using EAL. 
Calculated natural frequencies are given in 
Table 1. Some of the correspondi ng calculated 
node lines are shown in Fig 11 

A flutter analysis software system, known 
as FAST 6 , was used to predict the aeroelastic 
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Fig. 9.- Finite element model of the 
experimental wing. 


instabilities of the model. FAST calculates 
planar subsonic kernel function lifting-surface 
aerodynamics based on the natural mode shapes, 
natural frequencies, and generalized masses (as 
calculated through EAL) for the wing being 
analyzed. Four primary natural vibration modes 
(identified in Table 1, excluding first 
chordwise bending) were used in the present 
analysis. An aeroelastic instability solution 
is then completed by FAST with these calculated 
aerodynamics by the k method 7 . Divergence 
calculations are made by FAST provided the 
aerodynamics are calculated for a reduced 
frequency of zero. 

Aeroelastic calculations were made for the 
various sweep angles by varying the sweep of the 
analytical lifting surface. In the experiment, 
the wing was pivoted to the various sweep angles 
about a point on the trailing edge of the wing 
at the wing root. Therefore, a small portion of 
the wing near the root constraint was not in the 




(b) Torsional stiffness 
Fig. 10.- Measured stiffnesses. 


flow (due to the splitter plate, see Fig. 4) In 
the forward-sweep positions. The planforms used 
in the analysis for the various sweep positions 
corresponded to the appropriate portion of the 
experimental model exposed to the flow. A small 
gap between the splitter plate and the wing 
surface allowed the structure to be identical in 
all the sweep angles Therefore, the simulated 
structure of the model was identical for all of 
the sweep configurations. In addition, the 
structure was considered to remain the same for 
the two wing tip shapes. The physical change in 
the wing tip shape had no significant affect on 
the structural properties The wing tips were 
modeled analytically by changing the planform 
geometry to account for the aerodynamic changes 
induced by the wing tip shapes. The analytical 
wing tip shapes correspond to the geometry shown 
in Fig. 5 without the airfoil -shaped body-of- 
revolution fairings on the extreme tip. 





(c) First torsion H) Third hendino 

f a = 62 ^ Hz f a = 9? 3 Hz 

Fio. 11.- Calculated node lines and natural 
frequencies. 


RESULTS AND DISCUSSION 
Subcritical Response Techniques 

Divergence predictions were made during the 
wind-tunnel test using the improved static 
Southwell method and the divergence index 

method. Typical basic data used in the 

divergence predictions is shown in Fig. 6 for 
the rectangular tip model. Fig. 7 and Fig. 8, 
respectively, show typical improved Southwell 
method and divergence index method divergence 
prediction measurements made during the wind- 
tunnel test. Although these methods generally 
provided linear data from which to make 
divergence predictions. Fig. 7b shows the type 
of experimental scatter to which such methods 
may be susceptible. In this particular case, 
A = -60°, the higher dynamic pressure data was 
given greater weight by ignoring the three 
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lowest dynamc pressure data points shown in 
Fig. 7b. The divergence index method, Fig 8b, 
was similarly affected. In this case, the first 
three data points also were ignored and the 
fourth dynamc pressure data point was used as 
the reference condition 

An indication of the possible accuracy of 
these two prediction methods is shown in Figure 
12a-e. These plots show how the predicted 
divergence dynamc pressure varied as the 

experimental dynamic pressure was increased. 
Both axes of the figure have been normalized by 
q re f which is the divergence dynamc pressure 
predicted by the Southwell method after data had 
been obtained at the highest experimental 
dynamic pressure. In order not to damage the 
model, the dynamic pressure was not increased 
until divergence actually occurred. Based on 
the reference condition chosen, Fig 12 shows 
that for predictions made at a dynamc pressure 
within 35% of the final predicted divergence 
condition, the predictions are accurate within 
about ±10%. On the other hand, Fig. 12 also 
indicates the lack of extreme accuracy because 
of the experimental scatter exhibited by this 
model . 


Test Results 

Divergence measurements were made during 
the TOT wind-tunnel test at M = 0.4 for sweep 
angles of 0°, -15°, -30°, -45°, and -60° with 
the rectangular tip shape A tip parallel to 
the flow was tested only in the A = -45° 
position. The experimental results obtained 
with the Southwell method are shown in Fig. 
13. These results show the dramatic decrease in 
divergence dynamic pressure that occurred even 
at small forward sweep angles compared to the 
unswept wing and that the most critical 

configuration was at A - -45°. 

The effect of the two tip shapes tested is 
also shown in Fiq. 13. The tip parallel to the 
flow increased the divergence dynamc pressure 
by 14 percent compared to the rectangular tip at 
A = -45°. Tip effects similar to these have 
been documented previously for much lower aspect 
ratio winos 8 . These two tips are identical in 
the unswept position, so the effect of the tip 
may be proqressi vely less as the sweep anqle is 
decreased from A = -45° to A =0°. 


Analytical and Experimental Results 

The FAST aeroelastic software system was 
used to calculate divergence conditions for the 
model in each sweep configuration Four natural 
modes were used in the analysis. These modes 
were the three lowest frequency vertical 
(flapwise) bending modes and the first torsion 
mode {Fig. 11). Analysis also was made for the 
parallel tip shape in the A = -45° 

configuration The tip shape was modeled in the 
analysis by adjusting the planform geometry so 
that the tip was either rectangular with respect 
to the wing planform or parallel with the 
free-stream flow direction. The correlation 
between experiment and analysis at M = 0.4 can 
been seen in Fig. 13. The analysis was 
significantly conservative at all forward sweep 
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Fiq. 12.- Variation of predicted diveroence 
dynamic pressure as diveraence is 
approached. 
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^iq. 13.- Divergence boundaries for both tip 
models compared with calculated 
boundaries. 



(a) A = 0°, M = 0.5, p = 0.00207 sluas/ft 3 



angles. At A = -45°, the analysis was 14 
percent conservative. In the unswept position, 
analysis matched the experimental results. The 
effect of the two tip shapes tested at A = -45° 
was predicted accurately by analysis. Typical 
calculated plots of velocity-versus-incremental 
damping are presented in Fig. 14a-b. Divergence 
was the critical instability at all sweep 
angles, although flutter was analytically 
predicted to occur at nearly the same dynamic 
pressure as divergence at A = 0° (Fig. 14a). 
This predicted flutter condition is primarily a 
coupling of the second vertical (flapwise) 
bending mode and the first torsion mode. 
Similar flutter results have been shown to occur 
in other forward-swept wing experiments (see 
ref. 4). In some instances, flutter is actually 
the critical instability at small forward sweep 
angles (-15° < A < 0°). This was not the case 
for the high-aspect ratio wing tested in this 
experiment. For this model, as the sweep angle 
was decreased, the flutter condition rapidly 
moved to higher velocities (higher dynamic 
pressures). Fig. 14b shows clearly that 
divergence is the predicted instability at 
A = -45°. 

Analysis was conducted at M = 0.3, 0.5, and 
0 6 to determine the Mach number effect on 
divergence. The analytical divergence results 
are shown in Fig. 15. These results indicate 
that the divergence dynamic pressure decreases 
with increasing Mach number. The slope of this 
trend is very small at A = -45°, but increases 
as the sweep angle is moved away from the 
A = -45° position. 



Fia 15.- fa 1 cul effect of M ac^ nunher on 
divergence for various swppd annl^s 


( K) A = -AF° 5 M = n 5, p = O.wnf? sluos/ft 3 

Fi a. 14 - Results of the aeroelastic analysis 
near matched point conditions 
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CONCLUDING REMARKS 

A study of static aeroelastic divergence 
has been conducted for a high-aspect ratio, 
forward-swept wing model. Divergence 

characteristics were determined experimentally 
and analytically at sweep angles of 0°, -15°, 
-30°, -45°, and -60°. A rectangular tip shape 
was tested for all sweep angles. In addition, 
for A = -45°, a tip parallel to the flow was 
tested. 

The experimental divergence conditions were 
predicted using subcritical response test 
techniques. The diverqence dynamic pressure 
predictions improved as the divergence condition 
was neared. 

The wing tip shape had a significant effect 
on the divergence phenomenon For A = -45°, a 
wing tip parallel to the flow increased the 
divergence dynamic pressure by 14 percent 
compared to a rectangular wing tip. 

The analytical predictions using subsonic 
lifting surface theory were significantly 
conservative. At A - -45°, the analysis 

predicted a divergence dynamic pressure 14 
percent lower than the experimental value. 
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